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The distribution of Io’s volcanic activity probably reflects the position and 
magnitude of internal tidal heating, but full mapping of the moon has not 


been possible due to a lack of polar coverage. Here we use new observations 
of lo's polar regions by the Juno spacecraft Jovian Infrared Auroral Mapper 
to complete near-infrared coverage, revealing the global distribution and 
magnitude of thermal emission from lo’s currently erupting volcanoes. We 
show that the distribution of volcanic heat flow from 266 active hot spots 

is consistent with the presence of a global magma ocean and/or shallow 
asthenospheric heating. We find that Io’s polar volcanoes are less energetic 
but approximately the same in number per unit area than those at lower 
latitudes. We also find that volcanic heat flow in the north polar cap is 
greater than that in the south. The low volcanic advection seen at lo’s poles is 
therefore at odds with measurements of background temperature, showing 
that lo’s poles are anomalously warm. We suggest that the differences in 
volcanic thermal emission from lo’s poles compared with those at lower 
latitudes are indicative of lithospheric dichotomies that inhibit volcanic 
advection towards lo’s poles, particularly in the south polar region. 


The extreme level of volcanic activity on lo, the most volcanically active 
object in the Solar System’, is the result of tidally induced internal heat- 
ing’. Models predict enhanced heat flow at the poles of lo if tidal heating 
is deep inthe mantle and at lower latitudes if heating is predominantly 
in the asthenosphere, or a magma ocean is present® ’. Although the 
volcanoes of lo have been observed, and thermal emission quantified, 
at infrared wavelengths for decades", global mapping of volcanic 
activity has not been possible until now. Telescopic observations 
and previous flybys by spacecraft equipped with short-wavelength 
infrared imagers were mostly confined to viewpoints in the equatorial 
plane' that afforded poor polar coverage. In contrast to these previous 
observations, when Juno entered into a polar orbit of Jupiter in 2016", 
imaging of the lo’s poles became possible by the Jovian Infrared Auroral 
Mapper (JIRAM) at 3.5 pm and 4.8 um at spatial resolutions (as of July 
2022) of 151 km per pixel to 20 km per pixel. 


Although numerous hot spots on Io have been previously identi- 
fied in JIRAM data”, we do not draw any conclusions regarding the 
heat flow of lo from the number of hot spots alone. As the integrated 
thermal emission from individual hot spots on lo spans more than six 
orders of magnitude’, estimates of heat flux are more representative 
of volcanic advection than the hot spot number ina given area. Arecent 
study has determined that 4.8 um spectral radiance is a reasonable 
proxy for hot spot total thermal emission”. Although that study used 
Galileo Near-Infrared Mapping Spectrometer (NIMS) data obtained 
from Io observations almost entirely in the equatorial plane (so with 
poor viewing of the poles), an examination of a selection of lo hot 
spots observed by JIRAM at both 3.5 um (L-band) and 4.8 um (M-band) 
shows that polar hot spots are no different in colour temperature 
from those at lower latitudes”. The implications of this are discussed 
subsequently. 
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Fig. 1| Hot spot detections. The maximum, unsaturated 4.8 um (M-band) 
spectral radiances from 266 hot spots identified in Juno JIRAM data obtained 
from March 2017 to July 2022, using data from orbits PJ5 to PJ43 grouped by order 
of magnitude. The larger the symbol, the greater the 4.8 um spectral radiance. 


>1-10 GW um” >10-100 GW um! 

This is an equal-area Mollweide projection centred on 180° W, 0° N. The grid size 
is 30°. 4.8 um is a wavelength sensitive to the thermal emission from young, hot 
lava on the surface of 10%. Our hot spot detection numbers and distributions 
(Tables 1 and 2) differ substantially from those of ref. 17. 


Here we report the quantification of 4.8 um spectral radiance from 
266 individual hot spots identified in the JIRAM data (Fig. 1) obtained 
between 27 March 2017 (orbit PJ5) and 5 July 2022 (orbit PJ43). This 
wavelength band is highly sensitive to thermal emission from young 
lava surfaces on lo for the eruption temperature (~1,430 K for basalt) 
down to =200 K (ref. 20) (=70 K above the diurnal peak background 
surface temperature). 


Results 

We use, in part, superposition of JIRAM images (Fig. 2) obtained during 
each orbital encounter to reveal the presence of low spectral radiance 
hot spots that otherwise would not be easily seen in single frames. The 
techniques used to reduce the data are described in Methods. Our 
set of detected hot spots overlaps only 64% of those reported from 
another analysis of the same data”, but we do find many others (see 
‘Other hot spot detections’ in Methods). This is possibly the result 
of our geolocation technique creating superposition products with 
greater positional accuracy and sensitivity to thermal emission from 
faint hot spots than in the previous study”. The smallest 4.8 um spec- 
tral radiance we detect is 0.0091 GW um” in PJ41 data. This radiance 
is from a patera designated PV32 (JRM183), located at 74.5° W, 34.5° N 
just south of Zal Patera. Hot spot locations and maximum unsaturated 
4.8 um spectral radiances, maximum 4.8 um spectral radiances and 
average 4.8 um spectral radiances are provided in Supplementary 
Table S1. Table 1 summarizes our results. We define polar regions as 
being at latitudes >60° because models of interior tidal heating that 
calculate radially integrated heat flow** predict enhanced endogenic 
heat flow above latitude 60° if tidal heating is deep-seated in the mantle 
of lo, and less heat flow at these high latitudes and higher heat flow at 
lower latitudes if tidal heating is mostly asthenospheric or ifa magma 
oceanis present’. We find that polar regions collectively have a slightly 
lower hot spot density per unit area (hot spots per 10° km’) than those 
at lower latitudes. However, 4.8 um spectral radiance density per unit 
area (GW pm km”) from the polar caps of Io is <50% of that seen at 
lower latitudes, even when excluding Loki Patera, a powerful thermal 
source 180 km in diameter located at 310° W, 12° N. Loki Patera is a 
unique feature on Io in terms of size, persistence and magnitude of 
thermal emission’™. Previous studies have often treated Loki Patera 
as an outlier, running two analyses—one including Loki Patera and one 


excluding Loki Patera”. Excluding Loki Patera makes no significant 
difference to our conclusions. 

Furthermore, there is a difference in hot spot density and spectral 
radiance density between the polar caps. We find fewer active volca- 
noes inthe south polar cap than in the north. This difference, however, 
may be the result of south polar observations being fewer and at lower 
spatial resolution than those at the north polar cap (Figs. 2 and 3). 
The number of hot spots is not as important as how much energy is 
being delivered to the surface. In terms of 4.8 um spectral radiance, 
the south polar volcanoes generate a spectral radiance per unit area 
of 7kW um? km”, only half that seen in the north (15 kW um? km?) 
and only a quarter of that seen at lower latitudes (26 kW um“ km”). 
Although higher-resolution observations of the south pole may reveal 
additional hot spots, these are likely to have small outputs as they have 
so far been undetected. If this is the case, then these hot spots will not 
significantly alter the south polar value of spectral radiance per unit 
area. On average, the polar volcanoes of lo individually generate less 
energy than volcanoes at lower latitudes; and the south polar volcanoes 
generate less energy per volcano than the north polar volcanoes. 

In addition to using maximum unsaturated spectral radiance val- 
ues (Table 1), we also performed the same calculations using saturated 
pixel values (Table 2 and Methods). The spectral radiance totals and 
spectral radiance densities for hot spots in which saturated values are 
used are therefore minimum values. Two south polar cap hot spots 
(Tiwas Patera-JRM191 and Upulevo Fluctus-JRM214) and two north 
polar cap hot spots (Tvashtar A-JRMO84 and Lei-Kung C-JRM185) have 
high values owing to saturation, but we find that many lower-latitude 
volcanoes exhibit similarly high saturated radiances, so the same 
conclusions regarding distribution of spectral radiance are drawn. 

In further sensitivity analyses, we reproduced the hot spot and 
spectral radiance densities using averaged unsaturated radiance val- 
ues; we examined the variability of the spectral radiance quantities 
as a function of where the polar cap is defined, calculating hot spot 
spectral radiance density values across a 5° wide band either side of the 
60° latitude boundary. We also ran a model that excluded detections 
in which the emission angle exceeded 75° to ensure that low-latitude 
volcanoes were not generating higher spectral radiances than polar cap 
volcanoes simply owing toa high-emission angle cosine correction. We 
find no effect on our conclusions. 
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Fig. 2|JIRAM 4.8 pm observations of Io. a,c, Highest-resolution JIRAM 
observation of lo's south polar region during orbit PJ10 (112 km per pixel) (a) 
and north polar region during orbit PJ43 (21 km per pixel) (c). b,d, The PJ10 data 
overlain on the Galileo/Voyager base map” with some prominent volcanoes 
identified (b) and the PJ43 data overlain on the same base map (d). Sub- 
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spacecraft points are 317° W, 78° S for PJ10 and 175° W, 79° N for PJ43. The red plus 
is the respective location of the north or south pole. Differences in appearances 
between a and c are due to both the 4.3-fold factor in spatial resolution and 
differences in exposure time leading to more detector saturation in PJ10 images. 


Table 1| Distribution of hot spots and 4.8 um spectral radiance (orbits PJ5 to PJ43) using maximum unsaturated values and 


using a 60° latitude cap 


Region of lo Latitude Area(km?) Number Hotspot Total4.8um Mean hot Standard Skewness Standard 4.8um 

range (°) of hot density spectral spot4.8um deviation in 4.8um errorin spectral 
spots (number radiance spectral 4.8 um spectral 4.8 um radiance 
per (GWum") radiance spectral radiance spectral density (kW 
10°km?) (GWum") radiance (GWum") radiance yum km?) 
(GWum") (GW um") 

North polar cap 60-90 2.7710 20 7.21 40.79 2.04 5.03 3.17 112 14.714+5.32 

South polar cap -60 to 2.77108 12 4.33 19.80 1.65 2.03 1.56 0.59 714+2.58 
-90 

Both polar caps 60-90 5.54x10° 32 5.78 60.59 1.89 412 4.22 0.33 10.94 +1.97 
and -60 
to -90 

Lower latitudes <60and 3.59x10" 234 6.52 946.43 4.04 8.08 4.43 0.53 26.39+0.74 
>-60 

Lower latitudes, <60and = 3.59x10 233 6.49 877.91 3.77 6.65 3.74 0.45 24.45 +0.68 

excluding Loki >-60 

Patera 

Global 90 to 414x10" 266 6.43 1,006.83 3.80 7.74 4.30 0.47 24.32+0.59 
-90 

Discussion 


We have created a global map of 4.8 um spectral radiance that has good 
sensitivity to high latitudes from 266 volcanic hot spots identified in 
Juno JIRAM” data. We find 60% more volcanic spectral radiance per unit 
area at latitudes lower than 60°. Although the numbers of hot spots per 
unit area are not greatly different, polar volcanoes emit less energy at 


4.8 um than volcanoes at lower latitudes by a factor of more than two. 
This result is consistent with models of a global magma ocean or tidal 
heating in the shallow asthenosphere. Our results also show that from 
the analysis of currently available JIRAM data, north polar volcano spec- 
tral radiance is more than twice that of south polar volcanoes, suggesting 
dichotomies in structure and volcanic advection between polar regions. 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-023-02123-5 


Table 2 | Distribution of hot spots and 4.8 um spectral radiance (orbits PJ5 to PJ43) including saturated values and using a 


60° latitude cap 
Region of lo Latitude Area(km?) Number Hotspot  Total4.8um Meanhot Standard Skewness Standard 4.8um spectral 
range (°) of hot density spectral spot4.8um deviation in 4.8um errorin radiance density 
spots (number radiance spectral 4.8um spectral 4.8um (kWum”km?) 
per (GWum") radiance spectral radiance spectral 
10°km’) (GWum") radiance (GWum") radiance 
(GWum") (GWum") 
North polar cap 60-90 2.77 10° 20 7.21 58.67 2.93 6.66 3.10 1.49 21.18+7.65 
South polar cap -60 to 2.77*10° 12 4.33 44.52 3.71 6.37 2.32 1.84 16.07+5.80 
-90 
Both polar caps 60-90 5.54x10° 32 5.78 103.19 3.22 6.46 2.43 1.02 18.63+3.36 
and -60 
to -90 
Lower latitudes <60 and 3.59x107 234 6.52 1,388.09 5.93 15.69 TAS 1.03 38.71+1.08 
>-60 
Lower latitudes, <60 and 3.59x10’ 233 6.49 1,233.96 5.30 11.79 7.35 0.81 34.37+0.96 
excluding Loki >-60 
Patera 
Global 90 to 414x10” 266 6.43 1,494.39 5.62 14.90 712 0.91 36.02+0.87 
-90 
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Fig. 3 | Surface coverage by JIRAM. JIRAM surface coverage from orbits PJ5 to PJ43. The coverage of the south polar region is less frequent than that of the north 


polar region. 


Our analysis reveals a global snapshot of the locations where the 
volcanoes oflo are in active eruption, where new lava is currently being 
emplaced onto the surface as lava flows or exposed in overturning lava 
lakes’. We do not see thermal emission from areas that were active 
within the past few decades and centuries where surface temperatures 
have dropped below the JIRAM detection limit, around 180 K. The total 
4.8 um spectral radiance from the 266 hot spots within this dataset is 
1.5 TW pm ‘using maximum measured radiance (including saturated 
pixels) and1 TW pm ‘using maximum unsaturated radiance. Using the 
empirical relationship between 4.8 um spectral radiance and inferred 
total hot spot thermal emission described earlier”, these numbers yield 
estimates of total thermal emission between 25 TW and 18 TW. This 
compares favourably witha previous analysis of 4.8 ym ground-based 
telescope data, which yielded a global volcanic thermal emission esti- 
mate of 21 TW (ref. 23). 

JIRAM so far has therefore detected =30 TW (=55%) of the esti- 
mated 56 TW emanating from all of the volcanic edifices of lo'***. The 
JIRAM-derived volcanic thermal emission estimate is =28% of the total 
global thermal emission of lo (106 TW)”. 

Volcanic 4.8 um spectral radiance has been shown to be (and used 
as) a proxy for total hot spot thermal emission from the analysis of 


Galileo NIMS data’*”’, and this study shows that polar volcanoes of lo 
are not particularly different from those at lower latitudes in all aspects 
except areal extent. If lo were being predominantly tidally heated in the 
deep mantle, then the expectation is that polar volcanoes of Io would 
reflect this by being more numerous, or larger and more energetic, 
and having higher temperature lavas”. Regarding the first term, we 
find slightly fewer hot spots in polar regions. Nor do we find that the 
detected hot spots are larger in terms of power output than those at 
lower latitudes to make up for any scarcity at higher latitudes. 

The two latter terms—energy and temperature—are harder to 
disentangle. Considering JIRAM imager wavelengths, it is important 
not to interpret a higher colour temperature as evidence of higher 
magma eruption temperature. JIRAM imager data cannot be used to 
tightly constrain lava eruption temperature (and hence lava composi- 
tion) from measurements of thermal emission. Such a constraint can 
only be imposed using a shorter wavelength (0.4-1.5 um) imager”. 
The shape of the near-infrared thermal emission spectrum is mostly 
determined by the manner in which lava is erupted; a lava fountain 
has a high colour temperature (which is itself an approximation of a 
complex surface temperature—area distribution), whereas a lava flow 
of the same composition lava with a well-developed surface crust 
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Fig. 4 | Extensive volcanic activity at Tvashtar Paterae. a, JIRAM super- 
resolution image coverage of Tvashtar Paterae (centre of image) and the 
surrounding region from PJ43 (5 July 2022). These images are centred on 124.3° W, 
62.3° N. At 20 km per pixel, details of the surface distribution of thermal emission 
at Tvashtar Paterae become apparent. b, The outlines (solid white line) of the 
prominent thermal sources overlain on the lo global mosaic”. Tvashtar Paterae is 
the site of episodic, vigorous volcanic activity, including lava fountaining feeding 
the emplacement of lava flows and possibly lava lake overturn?®”. 


has a much lower colour temperature. Nevertheless, we do not find 
evidence of higher colour temperatures at polar-region volcanoes, 
suggesting that eruption styles are not greatly dissimilar to those at 
lower latitudes. It is likely that, before the Juno encounters, analy- 
ses examining high-latitude eruptions were skewed by poor spatial 
resolution and restricted viewing angles, and where detections were 
subsequently dominated by vigorous eruption styles, such as lava 
fountaining, which were highly energetic and which generated higher 
than average colour temperatures, while smaller, more numerous 
eruptions were not detected. 

An open question, however, is whether the distribution and 
magnitude of volcanic thermal emission itself is a proxy for global 
heat flow”. The polar temperatures of lo appear to be warmer than 
accounted for by solar insolation” and it is not known whether this 
results from high endogenic polar heat or from some property of the 
polar surface material with, for example, a high thermal inertia’. 
If high endogenic heat flow is the cause, it is also possible that this 
could promote more polar volcanism. Yet the polar volcanoes of lo 
are less powerful and slightly fewer than those at lower latitudes. 
It may be that this apparent conundrum will not be resolved until a 
spacecraft maps the global background endogenic heat flow pattern 
of lo, separate from heat flow at ~2% of the surface of lo occupied by 
active volcanoes”””°. 

The principal means of testing models of the interior heating 
and volcanic advection of lo comes from observing the distribution 
and magnitude of volcanic activity, thermal emission and lava com- 
position”. The coupling of lo to Europa and Ganymede via the orbital 
resonance’ has an important role in how the strength of the tidal dis- 
sipation (and by implication, interior state) has evolved. If lo oscillates 
into and out of the orbital resonance”, this has strong implications for 
the volcanism of lo. When the orbital eccentricity of lois damped, tidal 
heating and deep-mantle heating decrease. A lack of enhanced polar 
volcanism suggests that solid-body dissipation in the deep mantle of 
lois not taking place”. In addition to mapping the background thermal 
emission of Io, constraining the interior state of lo will probably depend 
on future spacecraft measurements of the strength of magnetic induc- 
tion, having separated magnetic induction of lo from effects caused 
by the lo plasma torus’, coupled to measurements of the liberation 
and diurnal tidal potential Love number of lo, k, (ref. 34). Excluding the 
deep-mantle heating model as a consequence of this analysis of JIRAM 
data means that a strong magnetic induction implies that a magma 


ocean is present, which would be supported by a high value of k, and 
large liberation amplitude™; whereas a weak magnetic induction and 
lowk, and liberation amplitude values imply solid-body dissipationin 
the asthenosphere””**», 

AsJuno obtains even closer observations of lo, the quality of ther- 
mal detections improves with increasing spatial resolution. We have 
seen how higher-resolution images of individual hot spots show more 
thermal structure. A good example is Tvashtar Paterae. JIRAM observa- 
tions obtained during PJ43 (SJuly 2022) show multiple hot spots within 
the two northern paterae of the Tvashtar Paterae complex (Fig. 4). For 
the purposes of quantifying thermal emission by volcano location, care 
is needed to determine whether these hot spots are separate eruptions 
or different parts of a single event. Another example is Prometheus, 
which appears asa single source in data up to PJ43, when two hot spots 
are resolved in higher-resolution images. 


Methods 

Data reduction 

JIRAM dataare processed, pixel by pixel, from band radiance (W m” sr”) 
to spectral radiance (W m° sr” um”) and corrected for distance from 
spacecraft to surface (that is, spacecraft altitude). Data are further cor- 
rected for emission angle, assuming Lambertian emission, to calculate 
surface-leaving radiance in W um” (reported in units of GW pm for 
convenience). The M-band (4.8 um) filter bandwidth is 0.4975 pm 
(ref. 16). The JIRAM instantaneous field of view Q = 2.37767 x 10% rad. 
For each pixel, therefore, the following steps are taken to determine the 
surface-leaving spectral radiance. An emissivity € of lis used. 


(1) Calculate pixel spatial resolution = R,patiai = OS, in which S is the 
spacecraft altitude (km). 

(2) Calculate pixel area Aix (km) = (QS)’. 

(3) Convert band radiance /,,,4, to spectral radiance spectrala 
(W m° sr” pm’) by dividing anaa by filter width. 

(4) Multiply /,pectra,a by pixel area in m? to convert to W sr pm”. 

(5) Multiply /.pectrara by TE to convert to W pm“, assuming Lamber- 
tian emission. 

(6) An emission angle e correction yields the surface-leaving spec- 
tral radiance in W pm™= I, pectra/COS(e). 

(7) The spectral radiance is summed for the pixels comprising the 
hot spot. 


Frame registration 

We initially investigated hot spot detection in image products navi- 
gated solely using NAIF SPICE kernels. We found that the products 
required additional positional correction. Using a superposition tech- 
nique allows identification of faint thermal sources that otherwise 
would be difficult to spot in individual JIRAM frames. 


(1) Using Juno mission NAIF SPICE kernels, lo is identified within 
each JIRAM image frame. 

(2) Positions are adjusted using limb-fitting and tying hot spots to 
known stable surface features (the former if lighting condi- 
tions allow). This is a vital step in identifying surface locations 
as SPICE kernels alone have insufficient precision. These 
adjustments improve the alignment of hot spots from image 
to image; reduce motion blur in the resulting summed images 
and make it easier to map hot spots to known surface features. 
For example, we incorporated 12 images into the PJ17 superpo- 
sition product. The range of adjustments needed to align all 
images was -1.25 pixels in the x direction and 0.75-2.75 pixels 
in the y direction. As the sub-spacecraft spatial resolution was 
=68 km per pixel (from a spacecraft altitude of 289,000 km), 
these adjustments are large distances on the surface of lo, 
translating to large differences in longitudinal position estima- 
tion at high latitude. 
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(3) For each on-moon pixel, pixel centre latitude, longitude, 
spacecraft altitude and emission angle are calculated. We cre- 
ate saturation masks for each JIRAM observation based on the 
detector sensitivity and spectral radiance for both L-band and 
M-band images so saturated pixels are flagged. 

(4) Each lo-containing frame in an orbit is reprojected to the same 
point perspective map projection at 10x spatial scale of the 
original data (5x for PJ41). 

(5) Asuperposition product is generated from the average of each 
of these frames for each orbit. 

(6) We examine each detection in each superposition image, being 
careful to avoid counting ‘detections’ caused by data artefacts 
(such as diffraction spikes). Latitude, longitude, size (number 
of pixels) and average band radiance are measured for each 
detection; using emission angle and altitude, spectral radiance 
(GW pm) is calculated for each hot spot. 

(7) For hot spots where radiance is detected in pixels and/or por- 
tions of pixels off the disk of lo, additional superposition images 
were generated with three-pixel offsets to the Xor Y adjust- 
ments used in step 2. This allows the total band radiance to be 
measured and the spectral radiance to be calculated for these 
high-emission angle hot spots. 


This methodology prevents misidentification of hot spots and 
results in a robust dataset. Detections are shown in Fig. 1 and listed in 
Supplementary Table S1. 


Background removal 

For hot spots identified in daylight, the reflected sunlight component 
is removed by subtracting the average radiance value of the pixels sur- 
rounding the hot spot from each hot pixel in the detection, a technique 
previously described and applied to NIMS data”. 


Hot spot detection 

No single threshold value of radiance was used to detect hot spots 
owing to differences in background noise and lighting condition, 
which could make detection more difficult when looking for fainter 
spots. In nightside/eclipse observations, a threshold brightest pixel 
50% greater than the background level was used. For hot spot detec- 
tions on the dayside, the threshold is lower, sometimes as low as 1.2x 
the background level. Additional criteria were applied to the faintest 
dayside detections. Was the spot seen as a hot spot at other times on 
the nightside or in eclipse? Is the proposed hot spot in an area of high 
albedo, such as Ra Patera and Acala Fluctus? If so, and other criteria 
are not met, then no detection is recorded. By adopting a conserva- 
tive approach, some very faint hot spot candidates are not included 
in the detections list. This approach does not adversely affect the 
conclusions of this paper as these hot spots are at the low end of the 
power spectrum and therefore contribute little to total 4.8 ym spectral 
radiance. 


Other hot spot detections 

We have compared our analysis with previous hot spot identifications 
by ref. 17 derived fromthe sameJIRAM dataset. Zambon et al. used hot 
spot count alone to conclude that the polar regions of lo were more 
volcanically active than at lower latitudes, the precise opposite of 
our findings. Possibly explaining this discrepancy, we find significant 
differences between the numbers and positions of our hot spot detec- 
tions andthe hot spot list of ref. 17. Over orbits PJ10 to PJ33, we identify 
240 hot spots. Of the 243 hot spots listed by ref. 17, we find 156 (Sup- 
plementary Fig. S1). Supplementary Fig. S2 shows the comparison of 
analyses of PJ17 data. We surmise that many hot spots detected by ref. 17 
are multiple detections of mostly high-latitude hot spots, and the 
superposition technique we use is more sensitive to fainter hot spots 
not detected by ref. 17. 


Uncertainty analysis 

JIRAM band radiance data are not provided with uncertainties. To 
calculate uncertainties in the spectral radiance density value, the 
mean, standard deviation and skewness of spectral radiance are 
calculated. These values are reported in Tables 1 and 2. Given that 
this non-negative set of data produces standard deviations greater 
than half of the mean values, this indicates that the data deviate sub- 
stantially from a bell-shaped (normal) curve. This is indeed the case, 
reflected in the large positive values of skewness. Rather than use 
the mean value and standard deviation to quantify uncertainty, the 
standard deviation and skewness are used to calculate the standard 
error, also reported in Tables 1 and 2. A value of twice the standard 
error is used to represent uncertainty in spectral radiance density in 
Tables 1 and 2. 


Saturated data 

JIRAM band radiance values plateau at the point of detector saturation. 
TheJIRAM detector saturates in a predictable manner. Supplementary 
Fig. S3 shows the saturation band radiance as a function of image inte- 
gration time. We identify and flag all saturated pixels and have created 
a saturation mask for every frame in every JIRAM observation of lo. 


Data availability 

The data analysed in this paper are available from NASA’s Planetary 
Data System at https://pds-atmospheres.nmsu.edu/data_and_ser- 
vices/atmospheres_data/JUNO/jiram_orbits.html. The products 
generated by this analysis will be available via the lo Geographical 
Information System Database” at Arizona State University (https:// 
rgcps.asu.edu/gis_data/) on publication. The lo Global Photomo- 
saic” used as a background image in Figs. 1, 2 and 4, and Supplemen- 
tary Fig. S3 is available from https://www.usgs.gov/media/images/ 
io-global-image-mosaic-and-geologic-map. 


Code availability 
Python scripts for data projection and for the generation of data back- 
planes, which for all pixels, ifthe pixel is on-planet, corrected latitude 
and longitude, band radiance, emission angle, spacecraft altitude, a 
saturation flag and spectral radiance, are available at https://github. 
com/volcanopele/juno. 
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